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[.  INTRODUCTION 


With  the  development  of  software  efficiently  implementing  quantum 
chemical  techniques  that  employ  energy  gradient  and  Hessian  data,  and  the 
concurrent  increase  of  access  to  modern  supercomputers,  systematic  studies  of 
molecular  structures  and  properties  have  become  routine.  Recently,  several 
groups  have  reported  efforts  to  compute  accurate  values  of  molecular  heats  of 
formation.  Pople  and  coworkers  described  a  study  of  first-row  hydrides  using 
full  fourth-order  perturbation  theory  with  a  variety  of  basis  sets.^  In  the 
same  proceedings,  Binkley  and  Frisch  described  results  of  perturbation  thepry 
studies  of  bond  dissociation  energies  for  a  series  of  first-row  diatomics. ^ 

The  results  described  in  these  papers  indicated  that  fourth-order  perturbation 
theory  calculations  using  extended  basis  sets  provided  generally  excellent 
results  for  sigma-bonded  systems,  with  less  accurate  results  obtained  for 
multiply  bonded  systems. 


Subsequently,  Page,  et  al . ,  used  a  similar  approach  to  evaluate  the 
dimerization  energy  of  borane  There  were  systematic  differences  among  these 
studies  that  should  be  noted.  Binkley  and  Frisch  computed  equilibrium 
structures  using  third-order  perturbation  theory,  whereas  the  other 
calculations  used  structural  parameters  predicted  at  the  Hartree-Fock  level. 
Page,  et  al.,  employed  an  extended  basis  set  when  computing  the  energy,  while 
the  other  authors  approximated  this  extended  basis  set  by  assuming  that  the 
effects  of  individual  basis  set  extensions  were  additive.  Related  to  these 
of  forts  to  predict  molecular  energies  accurately,  Handy  and  coworkers  have 
published  a  series  of  papers  that  attempt  to  evaluate  the  success  of  various 
levels  of  theory  in  predicting  the  properties  of  molecules Those 
computations  all  employed  extended  basis  sets  in  order  to  isolate  the  effect 
of  theoretical  method  from  the  effect  of  basis  set  deficiencies.  This  note 
describes  the  effects  on  computed  atomization  energies  of  variation  in  the 
theoretical  model  used  to  analyze  extended  basis  set  results.  In  particular, 
we  focus  on  differences  obtained  when  one  uses  either  the  isogyric  analysis 


employed  by  Pople 
previous  reports. 


et  al 

V  i* 


. ,  or  the  direct  method  of  computation  vised  in  our 
also  describe  the  application  of  isogyric  analyses  at 


several  orders  of  perturbation  theory. 


n.  MKTIIODS 


We  report  a  series  of  Mol  1 ei — Plesset  perturbation  theory  calculations  on 
small  molecules.  These  calculations  employ  1 arge  atomic  centered  basis  sets 
for  the  expansion  of  the  molecular  orbitals,  beginning  with  the  6-31 1G  basis 
designed  for  use  with  correlation  theory  methods. R  This  basis  is  augmented  by 
including  diffuse'  functions,  s-type  for  hydrogen  and  both  s-  and  p-type  for 
the  first  row  atoms.  Additionally  lhre<‘  si>ts  of  functions  in  the  first 
polarization  space*  and  one  seit  in  t  lieu  seu-einel  po  1  a  r  i  z.at  ion  space  are  aeided  for 
each  atom.  For  the  largest  me>lerule  studied,  diborane,  the  tevtal  number  of 
basis  functions  is  186.  Perturbation  theory  cal  ceilat  i  ons  wore*  performed  using 
the  GAUSSIANR2  computer  programs^  for  structures  opt.imizeei  at  the  SCF  level 
using  the  we’  1  I  -  kneie^n  demble  ze’ta  plus  polarization  basis  set.'^’*'  Structures 
for  those  rompoeinds  e:har;icter  i  zed  by  open-shell  wave’  func  t  i  ons  we’re  computed 
using  r  us  L  r  io  Led  Hartree-Fock  wave  func  t  ions  .  Harmonic  freque’neies  were  also 
computed  using  SCF  methods  with  the  I) 7P  basis. 


Results  are  reported  for  all  three  Bj  hydrides  and  diborane,  as  well  as 
water,  ammonia,  hydroxyl  radical,  carbon  monoxide,  boron  oxide,  BO,  and 
methane.  The  atomization  energies  of  the  compounds  are  evaluated  in  two 
ways.  One  technique  computes  directly  the  difference  between  the  sum  of  the 
energies  of  the  constituent  atoms  and  the  molecule’s  energy.  That  is,  we 
compute  the  energy  for  the  process, 

AH  >  A  +  nil 
n 

Alternately,  we  evaluate  the  energy  using  the  isogyric  technique 
descrihed  by  Pople,  et  al.,  in  Reference  1.  Tn  this  case,  energies  relative 
to  that  of  the  hydrogen  molecule  are  computed  using  reactions  for  which  the 
number  of  unpaired  spins  is  conserved.  This  is  equivalent  to  conserving  the 
number  of  electron  pairs  among  the  species  involved  in  the  comparison. 

Knowing  the  number  of  unpaired  electrons  for  a  particular  atomic  ground  state 
enables  one  to  write  isogyric  comparisons  with  molecular  hydrogen.  For  the  BH 
molecule,  one  has, 

BH  +  H  >  B  +  H2 

Assuming  that  theory  can  predict  the  energy  of  a  reaction  producing  an  atom 
plus  hydrogen  molecule(s),  then  the  precisely  Known  value  of  the  dissociation 
energy  for  the  ~  can  be  used  to  give  the  value  of  !>e  for  the  molecule  of 
interest. 

Our  computations  differ  from  those  of  Pople,  et  al  . ,  in  several  ways. 

The  theoretical  energies  used  in  their  analysis  were  obtained  by  improving  the 
results  of  energies  computed  using  fourth-order  Mol Ler-Plesset  theory  with  the 
6-11 1G**  basis  set.  Initially,  they  estimated  a  correction  to  the  correlation 
energy  by  using  an  extrapolation  method  that  estimates  the  contribution  made 
by  orders  of  perturbation  greater  than  four.  The  extrapolation  procedure 
produced  atomization  energies  that  differed  little  from  the  fourth-order 
results.  Tn  some  cases,  agreement  with  experiment  was  poorer,  an  indication 
that  the  main  source  of  error  was  the  incompleteness  of  the  basis  set  used  to 
desrihe  the  molecular  orbitals.  Additionally,  the  fourth-order  perturbation 
theory  calculations  Ignored  the  energy  contribution  due  to  excitations  of  the 
electrons  in  the  Is  core  of  the  first  row  atoms.  Most  of  our  calculations  do 
not  make  this  approximation,  but  neither  do  we  employ  the  extrapolation 
procedure  to  estimate  the  energy  contribution  obtained,  perhaps,  at  higher 
orders  of  perturbation  theory. 

The  basis  set  improvements  in  Reference  1  were  evaluated  incremental ly; 
the  effects  of  adding  diffuse  functions,  polarization  functions,  anil  functions 
in  the  second  polarization  space  were  assumed  to  he  additive.  In  the  present 
case,  all  calculations  are  done  with  the  extended  basis  set.  Pople,  et  al., 
concluded  that  the  addition  of  diffuse  functions  and  second  polarization 
functions  had  marked  effects  on  the  atomization  energy  calculation,  but  the 
effects  of  the  First  polarization  functions  were  primarily  equal  far  the  atoms 
and  the  molecules.'  If  the  effects  of  the  basis  set  extensions  were  truly 
orthogonal,  those  results  should  closely  mimic  the  data  descrihed  here.  Since 
it  is  unlikely  that  there  is  no  cooperative  effect,  our  approach  should  lead 
to  slight  differences  with  the  incremental  results. 


f. 


In  addition  to  the  full  fourth-order  perturbation  theory  calculations, 
the  diborane  calculations  were  performed  using  the  frozen  core  approximation. 
In  order  to  obtain  an  indication  of  the  effect  of  this  approximation  on  the 
computed  atomization  energy,  we  have  computed  energies  for  several  species 
using  both  a  full  perturbation  theory  calculation  and  the  frozen  core 
approximation.  In  addition,  we  have  included  the  third-order  perturbation 
theory  results  for  each  species.  We  anticipate  the  application  of  lower  order 
perturbation  theory  results  for  larger  molecules  using  semiempirical 
correction  factors.  Our  goal  here  is  to  provide  data  on  the  application  of 
this  more  approximate  level  of  theory  using  extended  basis  sets. 

III.  RESULTS 

We  summarize  the  results  of  the  perturbation  theory  calculations  in  Table 
l.  Our  experience  with  the  boranes  indicated  that  extended  basis  set 
calculations  yielded  atomization  energy  predictions  at  third-order  that  were 
nearly  equal  to  the  full  fourth-order  perturbation  theory  predictions.  With 
the  data  presented  in  this  table,  w.  can  compare  the  third-order  and  fourth- 
order  results  over  a  larger  class  of  compounds. 


Table  1 .  Elec 

■tronic  Energy 

Results  (Hartrees) 

lecule 

Hart  ree- Fock 

MI5  3 

SDT0-MP4 

SDT0-MP4(  F 

H 

-0  .49981 

-0  .49981 

-0  .4998  1 

B 

-24  .S3 108 

-24  .60644 

-24  .61167 

-24  .59308 

C 

-37  .69025 

-37  .79  163 

-37  .79638 

-37  .77642 

N 

-54  .39889 

-54  .52976 

-54  .53302 

-54  .51173 

0 

-74  .80934 

-74  .98919 

-74  .99316 

-74  .97094 

H2 

-1  .1  3275 

-1  .17027 

-1  .17178 

BH 

-25  .1  29  30 

-25  .23876 

-25  .24484 

8H~ 

-25  .75976 

-25  .87763 

-25  .88193 

bh3 

-26  .39969 

-26  .55272 

-26  .55784 

-26  .53676 

ch4 

-40  .21259 

-40  .4  5285 

-40  .46095 

-40  .43767 

nh3 

-56  .21978 

-56  .49086 

-56  .50067 

-56  .47654 

h9o 

-76  .05931 

-76  .35067 

-76  .3621  5 

-76  .3381  1 

OH 

-7  5  .419  82 

-75  .65542 

-75  .6621  5 

BO 

-99  .55687 

-99  .89308 

-99  .91786 

-149  .67687 

-1  50 .14.86  5 

-150.17612 

B'2H6 

-52  .83425 

-53  .1 3056 

-53.14242 

Table  2  contains  the  data  on  all  ten  molecules  studied.  Both  direct  and 
isogyric  predictions  are  reported  for  each  molecule?.  In  every  case,  the 
isogyric  computation  of  the  atomization  energy,  scaling  with  the  known 
atomization  energy  of  the  hydrogen  molecule,  predicts  a  larger  atomization 
energy  than  does  the  direct  method  of  computation.  The  difference  between  the 
direct  and  isogyric  comparison  is  equal  to  a  multiple  of  the  difference 
between  the  calculated  atomization  energy  of  the  hvdrogen  molecule  and  the 
known  value  that  we  employ  in  our  scaling. 
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Data  comparing  the  results  of  full  fourth-order  atomization  energy 
predictions  with  those  obtained  using  third-order  perturbation  theory  are 
presented  in  Table  3.  Both  direct  and  isogyric  results  are  considered.  The 
MP3  results  preserve  the  energy  difference  between  ttie  direct  and  isogyric 
methods,  although  the  magnitudes  are  greater  than  those  obtained  using  the 
fourth-order  data.  For  the  hydrogen  containing  molecules,  the  energy 
differences  between  the  third-  and  fourth-order  results  are  no1,  dramatically 
large,  nor  is  the  direction  of  the  energy  change  constant.  For  the  BO  and  O^ 
molecules,  on  the  other  hand,  the  energy  differences  are  large.  Since  neither 
molecule's  atomization  reaction  can  he  described  in  terms  of  the  breaking  of  a 
series  of  sigma  bonds,  it  is  not  surprising  that  the  reduced  treatment  of 
electron  correlation  provided  at  the  MP3  level  leads  to  so  different  a  result 
than  the  full  MP4  calculation.  Triple-excitation  diagrams  that  occur  first  at 
fourth-order,  are  required  to  describe  the  dissociation  of  these  multiply- 
bonded  systems  . 


Table  2.  MP4  Predictions  of  Atomization  Energies  (Units:  Mar  trees) 


Molecule 

Full 

Direct 

MP4 

Isogyric 

Frozen 

Di  rec  t 

Core  MP4 

Isogyr i 

BH 

0 .13336 

0.13666 

— 

bh2 

0 .27064 

0  .27294 

ch4 

0 .44674 

0  .45134 

0  .444  2  5 

0  .44884 

0 .66233 

0  .66923 

0  .66201 

0  .66891 

NH-j 

H~0 

0  .46822 

0  .4  7312 

0  .46  5  38 

0  .47228 

0  .36937 

0 .373«7 

0  .36755 

0  .3721  5 

oh 

0  .16918 

0.17148 

BO 

0 .31303 

0.31523 

— 

07 

0 .18980 

0.19210 

BA 

0  .95740 

0.96660 

Table  3. 

Comparison 

of  MP3  and  MP4 

Atomization  F.nergies  (Units:  H) 

SOTO- MP4 

MP3 

Mo  1  e nil  c 

Di rect 

T  sogy ric 

Di  rec t. 

Isogyr i r 

HI! 

0 .1 3336 

0.1  3566 

O.I325I 

0.1  3632 

BH., 

0  .27064 

0  .27294 

(I  .27157 

0  .2  75  38 

BH^ 

ch4 

0  .44674 

0.45  134 

0  .4468  5 

0  .43447 

0  .66233 

0  .66923 

0  .66198 

0  .67  3  4 1 

tl9o 

0  .46822 

0  .47512 

0  .46167 

0  .47310 

0  .369  37 

0  .37  397 

0  .36186 

0 .36948 

oil 

0 .16918 

0.17148 

0.16642 

0.17023 

BO 

0  .31303 

0  .31533 

i)  .29  745 

0 .30126 

n9 

0 .1  8980 

0 .19210 

0 . 1  702  7 

0 .17408 

H 


IV.  DISCUSSION 


Handy  and  coworkers  In  a  series  of  publications  investigated  the  effect 
of  the  level  of  correlation  energy  on  the  prediction  of  molecular  structures 
and  spectroscopic  properties  J*  1  Most  of  these  studies  used  an  extended 
basis  set  similar  to  that  employed  in  this  work.  In  each  case,  the  rationale 
is  that  the  extended  basis  set  removes  at  least  one  variable  from  the  analysis 
of  the  results;  additional  basis  set  improvements  will  have  scant  effect.  We 
assume  that  this  is  so  in  our  calculations.  While  we  have  not  set  out  to 
determine  the  convergence  of  the  perturbation  sequence  with  respect  to  the 
atomization  energy  predictions,  we  can  comment  on  that  based  upon  the  third- 
order  perturbation  theory  results. 

The  most  obvious  characteristic  of  these  calculations  is  the  systematic 
difference  in  predicted  atomization  energy  that  occurs  at  each  level  of  theory 
investigated.  Tsogyric  analysis  always  predicts  a  larger  atomization  energy 
than  does  a  direct  computation,  with  the  difference  between  the  methods 
increasing  as  the  completeness  of  the  theoretical  treatment  of  correlation  is 
reduced.  This  must  occur,  since  the  magnitude  of  the  correction  factor 
increases  as  the  level  of  sophistication  of  correlation  energy  calculation  is 
reduced . 

It  is  also  true  that  the  full  fourth-order  perturbation  theory 
atomization  energies  always  exceed  the  frozen  core  results.  The  differences 
are  generally  small,  a  consequence  of  the  fact  that  all  the  cases  for  which  we 
can  compare  numbers  involve  a  single  first-row  atom.  In  addition,  our  frozen 
core  predictions  are  always  equal  to  or  greater  than  the  predictions  given  by 
Pople,  et  al.,  for  the  same  molecules.  These  differences  are  primarily  due  to 
the  difference  in  hydrogen  atom  hasis  sets  used  in  the  two  sets  of 
ca 1 cul at  ions  . 

As  discussed  previously,  the  difference  between  the  direct  calculation  of 
the  atomization  energy  and  the  isogyric  prediction  is  a  multiple  of  the 
difference  between  the  computed  energy  for  the  hydrogen  molecule  and  the  exact 
value.  Thus,  the  correction  factor  for  the  BH  molecule  is  one-third  that  of 
the  methane  molecule.  To  evaluate  the  results  of  the  calculations,  we  compare 
the  atomization  energy  predictions  to  those  given  in  Reference  I .  Those 
authors  evaluated  zero-point  energy  contributions  using  the  best  vibrational 
frequency  data  available  at  that  rime?.  There  are  cases,  especially  for  the 
BH,,  molecule,  the  experimental  numbers  for  both  frequency  and  heat  of 
formation  are  clearly  incorrect.  These  empirically  derived  data,  along  with 
the  results  of  our  calculations  and  the  frozen  core  data  of  Topic,  et  al  .,  are 
summarized  in  Table  4.  Note  that  the  full  MP4  ca  1  cul  at  ions  using  the  isogyric 
method  of  analysis  overestimates  the  empirical  atomization  energy  in  every 
case  but  one,  the  BH,  molecule.  Our  value  lor  BH,  is  higher  than  that  derived 
i  n  a  recent  (experimental  study.  That  report  assigns  a  1)^  value  of 
26r>  .3  kcal/mole,  from  which  we  derive  a  D  equal  to  28H  ,U  kral/molo,  using  our 
previously  reported  zero  point:  energy.  Hur  direct  computation  of  the 
atomization  energies  tend  to  b<>  close  to  t  he  frozen  core  results,  and  art' 
always  lower  than  the  experimental  values,  except  for  the  BH  molecule.  Tn 
this  case,  there  is  reason  to  question  the  accuracy  of  the  theoretical 
result.  Both  the  BH  molecule  and  the  boron  .Horn  have  low-lying  excited  states 
that  are  the  consequence  of  double  excitations.  The  perturbation  theory 
techniques  used  in  this  study  do  not  accurately  account  for  the  energy 


contribution  due  to  these  excited  states.  Curtis  and  Pople  have  described  a 
technique  for  computing  corrections  due  to  higher-order  perturbation  theory 
terms, ^  and  they  have  demonstrated  the  application  of  the  technique  in  a 
study  of  small  boron  compounds.*^1  That  analysis  leads  to  a  reduction  by 
several  kcal/mole  of  the  atomization  energies  originally  reported  in  Reference 
1  . 

Since  the  magnitudes  of  the  energy  differences  between  theory  and 
experiment  are  so  small  it  is  not  possible  to  assign  the  "correct"  value, 
except  in  the  case  of  BH2  for  which  the  empirical  value  is  known  to  be 
incor rec  t . 


Table  4.  Theoretical  and  Experimental  Atomization  Energies  (kcal/mole) 


SDTQ-MP4 

1 

Molecule 

Direct 

Isogyric 

MP3 

SOTO- MP4 -FC1 

Experiment 

BH 

83  .7 

85.1 

85.5 

85  .3 

82  .8 

BH? 

169  .8 

171  .2 

172  .8 

168  .5 

196  .2 

BH“ 

280  .3 

283  .2 

285  .2 

280  .9 

280  .6  14 

ch4 

415  .6 

419  .8 

422  .5 

419  .9 

419  .8 

nh3 

29  3  .8 

298  .1 

296  .9 

296  .2 

297  .3 

h9o 

231  .8 

234  .7 

231  .8 

230  .9 

232  .2 

OH 

106  .2 

107  .6 

106  .9 

105  .5 

106  .6 

Finally,  we  consider  results  of  the  isogyric  analysis  with  the  third- 
order  perturbation  theory  results.  As  noted  above,  this  technique  fails 
utterly  when  applied  to  multiply  bonded  systems.  In  the  case  of  the  hydrides, 
however,  the  large  basis  set  predictions  are  in  reasonable  agreement  with  the 
full  fourth-order  results.  In  every  case,  the  corrected  third  order  results 
exceed  the  atomization  energy  predictions  of  the  direct  calculations  and  the 
frozen-core  calculations.  There  is  not  a  systematic  difference  between  the 
third-and  fourth-order  isogyric  results,  however.  While  admitting  that  the 
test  set  is  small,  we  must  conclude  that  the  third-order  perturbation  theory 
calculations  do  provide  surprisingly  reasonable  values  for  the  sigma  bonded 
molecules.  Certainly,  a  more  thorough  evaluation  of  the  application  of 
extended  basis  set  third-order  results  for  sigma  bonded  systems  seems 
warranted,  since  the  computational  advantages  of  this  lower  level  of  theory 
are  substantial.  The  recent  description  of  Carter  and  Coddard  of  an  excellent 
prediction  of  the  methylene  singlet-triplet  splitting  obtained  using  an 
extended  basis  set  with  the  computationally  efficient  CVB  method* ^  also 
indicates  the  potential  benefits  of  this  genre  of  quantum  chemical  technique. 
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